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I. INTRODUCTION 

In the Standard Model (SM), CP violation arises from 
an irreducible complex phase in the 3x3 quark mix- 
ing matrix V known as the Cabibbo-Kobayashi-Maska- 
wa (CKM) matrix [1, 2 . Unitarity of the CKM matrix 
requires that the relation V ud V* b + V cd V* b + V td V t * b = 0, 
which defines the unitarity triangle, be satisfied. The aim 
of the B Factories is to test the unitarity of the CKM ma- 
trix by the precise measurement of the angles and sides 
of the above triangle, whose nonvanishing area indicates 
violation of CP symmetry. 
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d > d d d 

(a) Tree (b) Penguin 

FIG. 1: Leading and sub-leading order Feynman graphs for 
the B° D c * )+ -D w ~decays. 

Both the BABAR and Belle collaborations have mea- 
sured the CP parameter sin2/3, where the angle f3 is de- 
fined as j3 = arg [— V c dV* b /V t dV t * b ]. The most accurate 
measurements of sin2/3 [HHS] use the b — > (ccjs transi- 
tion, in which B°'s decay to charmonium final states. 
Measurement of b —> ccd transitions such as B° —> 
£»(*)+£»(*)- s hould yield the same value of sin2/3 to the 
extent that the contributions from penguin processes may 
be neglected. 

The leading and sub-leading order Feynman diagrams 
contributing to B° — > L)(*)+£)(*)~ decays are shown in 
Fig. [l] The effect of neglecting the penguin amplitude 
has been estimated in models based on factorization and 
heavy quark symmetry, and the corrections are found to 
be a few percent [SI [7] . Loops involving non-SM particles 
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(for example, charged Higgs or SUSY particles) could 
increase the contribution from penguin diagrams and in- 
troduce additional phases. 

In T(4S) — s- B a B° events the time-dependent decay 
rate for B° — > D* + D*~ is given by 



P^(At) = 



[1 + Stag S n sin(Am ci Ai) 



4Tb l- ■ — 
5 t agCcos(Am d At)] , 



(1) 



where t\, is the B° lifetime averaged over the two mass 
eigenstates, Amj is the B°B° mixing frequency, and At 
is the time interval between the B° —> D* + D*~ decay 
(B lec ) and the decay of the other B (Btag) hi the event. 
The parameter Stag = +1 (—1) in Eq. (fTj) indicates the 
flavor of the Btag as a B° (B°), while 77 = ±1 indicates 
the CP eigenvalue of the B° D*+D*- final state. The 
parameters C and are given by 



C = 



1 + IAP 



23m(A) 



A = 



1± 
p A' 



(2) 



where A (A) is the matrix element of the B° (B°) decay 
and p and q are the coefficients appearing in the expres- 
sion of the physical mass eigenstates B^, Bh in terms of 
the flavour eigenstates B, B: 



\Bl) 
\Bh) 



: P \B) 

p\B) 



■q\B) 
q\B). 



Since the B° —> D* + D*~ is the decay of a scalar to two 
vector mesons, the final state is a mixture of CP eigen- 
states. The CP-odd and CP-even fractions have been 
previously measured from the angular analysis of com- 
pletely reconstructed events [HI [H] . 

A large deviation of the measured parameter in 
Eq.|2]fr om the value of sin2/3 measured in b — > (cc)s tran- 
sitions or a non-zero value of direct CP violation [T0HT2"] 
would be strong evidence of new physics. 

Both the BABAR [8 and Belle [5] collaborations have 
studied the CP asymmetries of B° — s- D* + D*~ decays 
using fully reconstructed events. In this article we report 
a new measurement based on the technique of partial 
reconstruction, which allows us to gain a factor of ~ 5 
in the number of selected signal events with respect to 
the most recent BABAR full reconstruction analysis in [S] . 
This result is complementary to the latter measurement, 
because the statistics used are largely independent of 
each other. 



II. THE BABAR DETECTOR AND DATASET 

The data sample used in this analysis has been col- 
lected with the BABAR detector [13] operating at the PEP- 
II asymmetric-energy B Factory located at the SLAC 
National Accelerator Laboratory. We have analyzed the 
full BABAR data set collected at the the T(AS) mass peak, 



\fs = 10.58 GeV, corresponding to an integrated luminos- 
ity of 429.0 fb" 1 . In addition, we have used 44.8 fb -1 
of data taken off-resonance to evaluate the background 
from events e + e~ — » qq, where q represents a u,d, s or 
c quark ( "continuum" ) . To study backgrounds and vali- 
date the analysis procedure, we use a GEANT4-based [14] 
Monte Carlo (MC) simulation in which coherent BB pro- 
duction is simulated using the package EvtGen [15j . 

The asymmetric energies of the PEP-II beams are an 
ideal environment to study time-dependent CP phenom- 
ena in the B°-B° system. The boost of the T(4S*) in 
the laboratory frame by f3j = 0.56 increases the separa- 
tion between the vertices of the two B meson daughters, 
allowing their precise measurement. 

The BABAR detector is described in detail in Ref. 13J. 
We give here only a brief description of the main com- 
ponents and their use in this analysis. Tracking is pro- 
vided by a five-layer silicon vertex detector (SVT) and 
a drift chamber (DCB). The SVT provides precise posi- 
tion measurements close to the interaction region that are 
used in vertex reconstruction and low-momentum track 
reconstruction. The DCB provides excellent momentum 
measurement of charged particles. 

Particle identification (PID) of kaons and pions is 
obtained from ionization losses in the SVT and DCB 
and from measurements of photons produced in a ring- 
imaging Cherenkov light detector (DIRC). A CsI(Tl) 
crystal-based electromagnetic calorimeter (EMC) enables 
measurement of photon energies and electron identifica- 
tion. These systems all operate inside a 1.5 T supercon- 
ducting solenoid, whose iron flux return is instrumented 
for muon detection, initially with resistive plate chambers 
and more recently with limited streamer tubes |16j . 



III. ANALYSIS METHOD 
A. Partial Reconstruction 

In the partial reconstruction of a B° — s- D* + D*~ can- 
didate, we reconstruct fully only one of the two D*^ 1 
mesons in the decay chain D* — > D°tt |17j . by iden- 
tifying D° candidates in one of four final states: Ktt, 
Ktttt , Kiririr, K^tttt. The vertexing algorithm fits the 
two-step decay tree simultaneously, correctly calculating 
correlations among all candidates. In the first three D° 
decay modes, assumed to represent Cabibbo-favored de- 
cays, charged kaon tracks are selected using PID infor- 
mation from the DIRC, SVT and DCB. In the last decay 
mode, candidates are selected by constraining pairs 
of oppositely charged tracks to a common vertex. 

Since the kinetic energy available in the decay D* — > 
D° 7r is small, we combine one reconstructed D* ± with 
an oppositely charged low-momentum pion, assumed to 
originate from the decay of the unreconstructed D* T , and 
evaluate the mass m ICC of the recoiling D° meson by using 
the momenta of the two particles. For signal events m rec 
peaks at the nominal D° mass [TH] with an r.m.s. width 
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of about 3 MeV/c 2 , while for background events no such 
peak is visible. Thus, m roc is the primary variable to 
discriminate signal from background. The calculation of 
77i rcc proceeds as follows (refer to Fig. [2] for definitions of 
the various momenta and angles that we use). 




FIG. 2: Momenta and angles in the T(4S) center of mass 
frame used in partial reconstruction. The orthogonal axes u 
and v are normal to the momentum po* of the reconstructed 
D* , and u lies in the plane defined by the momenta of the D* 
and slow pion, pr>* and p s . The angle <j> is measured in the 
u — v plane. 

The cosine of the angle between the momenta in the 
T(4S) center of mass (CM) frame of the B and the re- 
constructed D* is readily computed: 



cos6 BD * = 



-M g0 + EqmEd' 
2p B \pD*\ 



(3) 



where all particle masses are set to their nominal val- 
ues |18j . Ejy* and pp* are the measured energy and mo- 
mentum of the reconstructed D* in the T(4S) CM frame, 
-Ecm/2 is the energy of each beam in the CM frame, 

and pb 



— Mg is the B-meson CM momen- 
tum. Events are required to be in the physical region 
\cos6 B D' \ < 1. 

Given cos Obd* and the measured momenta of the D* 
and oppositely charged low- momentum pion, 7r s , the B 
four-momentum can be calculated up to an unknown 
azimuthal angle 4> around pw ■ For any chosen value 



of 4>, conservation laws determine the unreconstructed 
D a four-momentum </£>(</>), and one can thus compute 
the corresponding ^-dependent invariant mass m(<f>) = 
\/|to(^)| 2 ' The value of is not constrained by kine- 
matics and may be chosen arbitrarily, to the extent that 
the shape of the resulting m(0) distribution may still be 
described by the type of functions used in our fits. We 
have chosen the value for which cos0 — 0.62, which is 
the median of the corresponding Monte Carlo distribu- 
tion for signal events obtained using generated momenta, 
and define the recoiling D° mass m ICC = m(cos = 0.62). 
We use the same convention to obtain the direction of 
the unreconstructed D° meson. 



B. Backgrounds and Event Selection 



Backgrounds to the B° 
the following: 



D* + D* process include 



• Combinatorial BB background, defined as decays 
other than B° — > D* + D*~, for which the m rcc dis- 
tribution is approximately flat. 

• Peaking BB background, defined as decays other 
than B a -> D* + D*~, in which the m roc distribu- 
tion peaks in the signal region. It will be shown 
later that the contribution from this background is 
negligible. 

• Background from non-66 events. 

Combinatorial BB background events are reduced by 
the following requirements: For the K®irir mode, we re- 
quire the invariant mass of the pion pair to be within 25 
MeV/c 2 of the K ° mass [TH] . The corresponding vertex 
must be separated by more than 3 mm from the beam 
axis. For the Kirir mode, ir° candidates are formed 
from pairs of photons detected in the EMC, with energies 
greater than 40 MeV, for which the invariant mass differs 
by less than 20 MeV/c 2 from the nominal 7T° mass [18] . 
The reconstructed D° mass must be equal to the nominal 
one |18j within 2 or 2.5 standard deviations, depending 
on the D° reconstruction mode. The momenta in the 
T(4S) CM frame of the reconstructed D* and tt s from 
the missing D must be, respectively, in the range 1.3- 
2.1GeV/c and smaller than 0.6GeV/c. The difference 
AM = \Md* — Mjjo — Mk\ must be equal to the nom- 
inal [T5] value within 1 or 1.5 MeV/ c 2 , according to the 
presence or absence of DCH hits in the pion track ap- 
pearing in the reconstructed decay D* — > D° n. The 
probability of the vertex fits must be greater than 10~ 2 , 
for both the D° and the D* reconstruction. 

The requirement on the D° vertex fit probability intro- 
duces a small but measurable bias toward lower values of 
the B lifetime. Due to partial reconstruction, the tracks 
used to make the D° vertex may originate from the same 
or different B mesons. In the latter case, since not all 
tracks are from the same point in space, the \ 2 °f the 
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vertex fit tends to be bigger. This effect worsens with in- 
creasing distance between the two B decay vertices, caus- 
ing vertices further apart to be rejected more frequently. 
We have verified this on signal Monte Carlo events, for 
which we have measured a lifetime lower than the gener- 
ated value. Consequently, for the signal At probability 
distribution functions (PDF's) we use the value of 77, fit- 
ted to signal Monte Carlo. 

In events passing this selection we find more than one 
candidate decay chain in about 25% of the cases, usually 
differing only in the slow pion n s , but sometimes in the 
components of the reconstructed D* . When this hap- 
pens, we choose one candidate chain, based respectively 
on the largest number of DCH hits in the tt s , or accord- 
ing to a x 2 based on the reconstructed D° mass and AM 
quantity above. For signal Monte Carlo, the probability 
for this candidate chain to be the correct one is 0.95. 

The main suppression of continuum background is ob- 
tained by requiring that the ratio R2 of the 2 nd to the th 
Fox- Wolfram moment [19] , computed using all charged 
particles and EMC clusters not matched to tracks, be 
less than 0.3. 



C. Fisher Discriminant 

To further reduce continuum background, we com- 
bine several event-shape variables into a Fisher discrimi- 
nant |20) F. Discriminating power originates from the ob- 
servation that qq events tend to be jet-like, whereas BB 
events have a more spherical energy distribution. Rather 
than applying requirements on F, we use the correspond- 
ing distribution in the fits described in Sec. |III E] 

Our Fisher discriminant is a linear combination of vari- 
ables chosen, according to Monte Carlo studies, to maxi- 
mize the separation between BB and continuum events. 
The first nine variables describe the energy flow inside 
nine concentric cones centered around the direction of the 
reconstructed D* ± . In addition, we use the momenta of 
the charged and the neutral particle closest to the cone 
axis, the polar angles in the CM of the reconstructed D* 
momentum and the thrust axis T for charged tracks in 
the -Btag vertex (see next paragraph), the angle between 
the reconstructed D* momentum and T, and the sum 
S = SiPi x P2(cos 0i) over the B tag charged tracks, in 
which pi is momentum, Pi is the 2 nd Legendre polyno- 
mial of argument cos 9i, and 9i is the angle between track 
i at the origin and T. 

D. Flavor Tagging and Decay Time Measurement 

For this analysis, two measurements are needed: the 
difference At between the proper decay times of the par- 
tially reconstructed B meson and the other B meson in 
the event, and the flavor of the latter. 

The flavor tagging algorithm is based on tracks iden- 
tified as electrons, muons or kaons. The electron and 



muon tags contribute equally to the total sample and, 
since these events are kinematically almost indistinguish- 
able and have very similar effective tagging efficiency, we 
treat them as one homogeneous "lepton" sample. 

The tagging tracks must be chosen among those not 
used in B rcc reconstruction and must originate from 
within 4 mm (3 cm) of the interaction point in the trans- 
verse (longitudinal) view. The momentum of the lepton 
candidates is required to be greater than 1.1 GeV/c in or- 
der to reject most leptons from charmed meson decays. 
If one or more lepton candidates are qualified, the tag 
flavor is assigned based on the charge of the lepton with 
the highest center-of-mass momentum. If two or more 
qualified kaons are present, the event is used only if the 
flavor is unambiguous. If both a lepton and a kaon tag 
are available, the lepton tag is used. 

The time difference At is calculated using At = 
Az/7/3c, where Az — z rcc — z tag is the difference between 
the z-coordinates of the partially reconstructed B Tec and 
Stag vertices and the boost parameters are calculated 
using the measured beam energies. The uncertainty cta* 
on At is calculated from the results of the z rcc and zt ag 
vertex fits. We require |At| < 20 ps and a At < 2.5 ps. 

We define the B Tec vertex as the decay point of the 
fully reconstructed D* ± . The tt s track from the other 
D* ± is not used, since it undergoes significant multiple 
Coulomb scattering and hence does not improve the z rcc 
measurement resolution. 

The -Btag vertex reconstruction depends on the tagging 
category. For kaon-tagged events, we obtain z t ag from a 
beam spot constrained vertex fit of all charged tracks in 
the event, excluding those from the B rcc meson, and ex- 
cluding also tracks within 1 rad of the unreconstructed 
D° momentum in the CM frame, which presumably orig- 
inate from the D° decay. We require the probability of 
this fit to be greater than 10 -2 . For lepton-tagged events, 
we use the lepton track parameters and errors, and the 
measured beam spot position and size in the plane per- 
pendicular to the beams (the x-y plane). We find the 
position of the point in space for which the sum of the 
X 2 contributions from the lepton track and the beam spot 
is minimum. The z coordinate of this point is taken as 

-^tag • 

The beam spot is measured on a run-by-run basis 
using 2-prong events (Bhabha and fi + fj,~), and has an 
r.m.s. size of approximately 120 /im in the horizontal di- 
mension (x), 5 /im in the vertical dimension (y), and 
8.5 mm along the beam direction (z). The average B me- 
son flight distance in the x-y plane is 30 /im. To account 
for the B flight distance in the beam spot constrained 
vertex fit, 30 /im are added in quadrature to the effective 
x and y sizes. 



E. Probability Distribution Functions 

We use two PDF's, P on for on-resonance, and P g for 
off-resonance data. The former depends on the variables 
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rn rcc , F, At, (TAt, Stag, and is given by the sum of the 
PDF's for the different event types described above: 

Pan = J BB [/sig-fsig + (1 — /sig)-Pcomb] + (1 — /bb)-^?^) 

where P s ; g , P CO mb, and P qq are respectively the PDF's 
for signal events, for combinatorial background from BB, 
and for continuum. Moreover, f BB - is the fraction of BB 
events in our sample, and / s j g is the fraction of signal 
events in BB events. The PDF for off-resonance data, 
P s, is reduced to just one component, P qq , as the off- 
peak sample contains only continuum events. 

According to Monte Carlo, the distributions of B a B° 
and B + B~ combinatorial background events are very 
similar and can be described well by the same PDF. 

We do not consider the fraction of BB events a free 
parameter, but fix it to f BB - = 1 — f qq , where f qq is the 
fraction of continuum events in the on-peak sample and 
is defined by 



AW- 



off-peak ^on- peak 



on— peak ^off— peak 



(5) 



where A's are the number of events left by our selection 
in the on- and off-peak samples and £'s are the integrated 
on- and off-peak luminosities. 

Each of the Pi (i = sig, comb, qq) can be expressed as 
the product of three one dimensional PDF's: 



Pi(m rec , F, At, a At , Stag) = 

A^TOrec) Fi(F) T[(At, a At , S tag ), 



(G) 



that are the probability distributions of the recoil D° 
mass Mi(m rcc ), the Fisher discriminant function J-'i(F), 
and the decay time difference function T((At, aAt, St ag )- 
This follows from extensive Monte Carlo studies showing 
that the correlations among these variables are negligible. 



and F(F) PDF's 



The m rcc distribution of all sample components can be 
well modelled in the lower region of the spectrum with a 
so called "Argus function" [2"T] : 



A(m ICC ) = m Icc ^l - (m rcc /m cp ) 2 • ^'"W^ (7) 

where m op is the kinematic endpoint (m rcc < m ep ) and 
c is a free parameter describing the slope. This func- 
tion alone, however, is not sufficient to account for the 
abrupt fall of the m loc spectrum near the kinematic end- 
point. For the signal sample we model this shoulder with 
an asymmetric error function with widths cr; and oy, ta- 
pered off at low m rcc by an exponential factor with decay 
constant a: 



£(m r 
[ e m ' 



c/a [l - erf(m rcc - to cp )/(v / 2ct;)], m rcc < m c 
3 mrcc/a [l - erf(m roc - m op ) /(\/2oy)], m ICC > m e , 



Thus, we describe the signal m roc distribution with a 
combination of three functions: a Gaussian G having av- 
erage vtlq and standard deviation oq for the well recon- 
structed peaking component; an Argus function, mainly 
for events in which the right D* is combined with a low- 
momentum pion from another decay chain; and the £ 
function: 



A4 sig (m ri 



/ s f g • A(m mc ) + 

(1 - /sfg) • [f° ■ G(m r 

(1 - f G ) ■ £(m ICC ) . 



(8) 



In Eq.[8] f^ g is the fraction of events described by the 
Argus component and / is the fraction of events in the 
Gaussian peak relative to the non- Argus component. 

For the background, both combinatorial and contin- 
uum, we set the fraction of the Gaussian component to 
zero, and model the distribution at the endpoint with a 
simple error function of width a. However, for the case 
of combinatorial background in kaon-tagged events, we 
find that two different Argus components (_4i and A2) 
are needed to correctly describe the entire reconstructed 
mass spectrum. We thus define two PDF's according to: 



M comh (m rcc ) = /, 



comb ' erf ( m 

erf 



+ (1 /comb) 
+ (1 ~ /comb) 



^ep; ^comb) ~t~ 

■Ai(m Iec ) + 

A2(m iec 



J comb 



M qq (m rcc ) 



feri 
J qq 

(1- 



erf(m rec - 

/|| f )--AlKec). 



>);<7 qq ) + 



(9) 



(10) 



The parameter m op represents simultaneously the two 
Argus endpoints and the error function inflection point. 

The Fisher discriminant PDF Ti is parameterized by 
two Gaussian functions for each event type i = (BB, qq), 
having standard deviations and af~, and common 
mean /if. 



T.(P\ fexpRF-^m^) 2 ] F<^ 



Since the Fisher variable is designed to discriminate be- 
tween qq and BB events, we expect the Fisher discrim- 
inant for signal events to be indistinguishable from that 
of BB combinatorial events. We have verified this expec- 
tation with Monte Carlo studies, and thus use the same 
Fisher discriminant to describe both event types. 



2. At PDF's 

The Ai-dependent part of the PDF is a convolution of 
the form 



T/(At, CTAt, Stag) = 

J dAt tluc Tj(A<truej Stag) 72* (At — Aitruej ^At)) 



(12) 
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where T is the distribution of At trnc , the true decay time 
difference, and 72. is a resolution function that parame- 
terizes detector resolution and systematic offsets in the 
measured positions of vertices. 

Taking into account the mistag probability and the 
effect of tags due to the unreconstructed D°, the Ai tr ue 
signal PDF in Eq.[l2]can be written as 

T sig = 4^ e-l A **™l^ • {l - 5 tag Aw(l - a) + 

+ Stag (1 - 2w) (1 - a) (13) 
■[CcQs(Am d A* true ) + S sin(Am d Ai truc )]|, 

where the time-dependent CP asymmetry parameters S 
and C are the object of the measurement discussed in 
the present article and a (see Sec. IV B I is the fraction 
of events in which the tagging track is from the unrecon- 
structed D°. We parameterize possible detector effects 
leading to a small difference between the mistag proba- 
bility of B° tags (oj + ) and that of B° tags by using 
the average mistag rate uj = (us + +lj~)/2 and the mistag 
rate difference Auj = lu + —lu~ as parameters of the PDF. 

Since the BB combinatorial background is dominated 
by non-CP final states, the CP asymmetry is expected to 
be negligible. However, we allow the PDF to accommo- 
date some contamination from CP final states. There- 
fore, we parameterize the BB background At truc distri- 
bution with a PDF similar to that for signal events given 
in Eq. 13 We also add a fraction fg of a ^-function, to 
allow for a zero-lifetime component: 

/co m b^(|Aitruc|)( 

.fi _ t s \ . !_ 

V J comb / 4r com 
1 — Stag AoJcombH 



comb — 



Stag Awf omb 

-|At truo |/T comb 



(14) 



St: 



C comb COs(Am<jAt t rue) + 



Scomb sin(Am d Ai t ruo) | 



The second term of the PDF is obtained from Eq.[T3| 
with to — a — 0, as these are not defined for background 
events. The C com b, Scomb parameters describe small fluc- 
tuations in the Ai tmc distribution of background events 
and possible CP event contamination, leading to a small 
effective CP violation value. 

The Au parameters, which for signal events is the dif- 
ference in the mistag probabilities for B° and B°, allow 
for differences in the number of events tagged as a B° or 
B° in the same background sample. We use this PDF to 
describe both the B°B° and B + B~ components. 

The PDF for the background due to continuum events 
is modeled with a simple exponential decay distribution 
plus a fraction fg of a <5-function: 



T„. 



f- 

Jqq 



1 



Au/- 

99 



(1-4) 



1 - S, 



tagAWgg 



<S(|Att ruo |) + 

e -|At t 



i 



(15) 



where the parameters Aw^ and Auj qq - allow for differ- 
ences in the number of events tagged as a B° or B° in 
this sample. 



3. Resolution Functions 

The functions T[ of the measured time difference At, 
to be used in the fits, are obtained by convolving the Ti 
PDF's of Eq.[l3j[T4j[T5] with the appropriate resolution 
function for events of type i (i=sig, comb, qq). 

The resolution functions are parameterized as the sum 
of three Gaussian functions: 



Ki(t r ,(TAt) = fj 1 Qi(t r ,(TAt) 

+ (i-f?-msr(u,<TAt) 



(16) 



where t r — At — At trU c is the residual of the At mea- 
surement, and Qf, Qf, and Q° are the "narrow", "wide", 
and "outlier" Gaussian functions. The narrow and wide 
Gaussian functions incorporate information from the At 
uncertainty a At, and account for systematic offsets in the 
estimation of a At and the At measurement. They have 
the form 



GHtr,<?At) 




(17) 



where the index k takes the values k = n,w for the nar- 
row and wide Gaussian funcions, and b\ and s\ are pa- 
rameters determined by fits. The outlier Gaussian func- 
tion, describing a small fraction of events with badly mea- 
sured At, has the form 



1 



exp 



2«) 2 



(18) 



In all fits, the values of b° and s° are fixed to ps and 
8 ps, respectively, and are later varied to evaluate sys- 
tematic uncertainties. 



F. Analysis Procedure 



After the event selection described in Sec. lIII Bl is com- 
plete, the rest of the analysis proceeds with a series of 
unbinned maximum-likelihood fits, performed simultane- 
ously on the on- and off-resonance data samples and inde- 
pendently for the lepton-tagged and kaon-tagged events. 
The procedure can be logically divided in the following 
three steps, which we shall discuss in detail in the follow- 
ing paragraphs: 

1. In the first step we determine the signal fraction 
/si g in Eq.Bland the shape of Ai(m rec ) and F{F) 
in Eq.[6]for the different classes of events (signal and 
backgrounds, kaon and lepton tagging categories). 
This is done by fitting data with the PDF 



P;(m rec , F) = Mi(m iec ) Fi(F), 



(19) 
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ignoring the time dependence; we refer to this step 
as the kinematic fit. 

2. In the second step we determine the tagging dilu- 
tion due to wrong tag assignments. 

3. In the last step we perform the time-dependent 
fit to the data. We fix all parameter values ob- 
tained in the previous steps and use the full PDF of 
Eq.[6]to determine the parameters of the resolution 
functions, T/(At, a At, Stag), and the CP asymme- 
try values C, S of the signal and of the BB combi- 
natorial background component. 

The fitting procedure has been validated using both full 
Monte Carlo and, where the requested number of events 
would be too large, the technique of "toy" Monte Carlo. 
In a toy Monte Carlo, events are described by a small 
number of variables which are generated according to our 
PDF's. 



IV. RESULTS 

Event selection yields the numbers of events listed in 
the top two rows of Table |T] The third and fourth rows 

TABLE I: Event selection yield. The first uncertainty shown 
is statistical, while the second uncertainty on the number of 
continuum events accounts for a 1% relative uncertainty on 
the on-peak and off-peak luminosities. 





# of events 






kaon tag 


lepton tag 


on-peak 


61179 


20855 


off-peak 


1025 


51 


continuum 


9814 ±307 ±196 4 


188 ± 68 ± 10 


BB 


51365 ± 364 


20367 ± 69 


N 3i& 


3843 ± 397 


1129 ±218 



show the number of continuum and BB events calcu- 
lated, using Eq.[5] from the number of off-peak events in 
the second row. The numbers of signal events in the last 
line of the table are calculated using the signal fractions 
obtained from the kinematic fit described in the next sec- 
tion. 



A. Kinematic Fit 

We begin by fitting the shape of our signal, 
■A/f S i g (m roc ), using a large sample of Monte Carlo signal 
events. The parameters most relevant to determine di- 
rectly the signal fraction in the data, and consequently 
our final result for S and C, will be released again in the 



final kinematic fit. They are (refer to Eq.|8]): the Gaus- 
sian fraction fa, mean value mo, and standard deviation 
cr G , and are shown in the last section of Tabic [TT] 

Next we fit the Fisher Tqs and recoil mass M.„q dis- 
tribution to the off-peak data sample. As the number of 
off-resonance events selected in the lepton tagged sample 
is too small to yield convergence, we set the lepton tag 
sample parameters to the corresponding values obtained 
from the fit to the kaon tag sample. Due to the small 
continuum fraction in the lepton sample, we judge that 
this does not introduce any significant systematic effect. 
The J- q q parameters are fixed in all subsequent fits. 

We initialize the parameters of the BB combinatorial 
background PDF directly from the data, using a sam- 
ple of events in which the contribution of signal events 
is much reduced. We obtain this sample by combining 
a D* with a pion of wrong sign charge (WS sample). 
We have verified, both on Monte Carlo and in the m rcc 
sideband for data, (1.836-1.856 GeV/c 2 ), that the shape 
of the M(m ICC ) distribution for combinatorial BB back- 
ground is well described by that of the WS data sample. 

To evaluate a possible contribution from a peaking 
component in the BB background events, we have al- 
lowed the Gaussian fraction fa in Eq.[8] to float in a fit 
to a sample of Monte Carlo background events; this frac- 
tion is found to be 0.000 ± 0.002, and is therefore set to 
zero. 

Finally we fit the on-peak data sample, leaving as free 
parameters the fraction / s j g of signal events in the BB 
component, some of the shape parameters of the contin- 
uum and BB combinatorial background *M C omt>, some of 
the signal parameters in A^ s ig, and the shape parameters 
of the Fisher discriminant J- ~ B g. Table ITT| summarizes the 
results and provides information about which parameters 
are released in the fit (statistical uncertainties given) and 
which ones are taken from previous fits (no uncertainty 
given) . 

The final results of the kinematic fits for the kaon and 
lepton tagged sample are shown in Fig. [3] and Fig. [4] 

B. Determination of mistag probabilities 

A common problem of analyses using the partial recon- 
struction technique is that a fraction of the tracks used in 
tagging may belong to the unreconstructed D , leading 
to a mistag of the event. As the tracks originating from 
the missing D° tend to align to its direction of flight, this 
fraction can be reduced by applying a constraint on the 
cosine of the CM opening angle 6* tag between the tagging 
track and the direction of the unreconstructed D°. We 
require cos# t ag < 0.75 (0.50) for the kaon (lepton) tagged 
sample. 

We find from signal Monte Carlo that before this re- 
quirement 26% (13%) of kaon (lepton) tags originate from 
a missing D°. We call these events "D-tags", while "T- 
tags" are those in which the tagging track (either direct 
or cascade) is from the tag B. 
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TABLE II: Results of the final kinematic fits. The values of fixed parameters are given without uncertainties. 



PDF 


parameter 


description 


kaon tags 


lepton tags 




/sig 


Signal fraction 


(7.5 ±0.7)% 


(5.5 ±1.1)% 




f^BB 




0.723 ±0.005 


0.721 ± 0.009 




BB 




0.361 ±0.003 


0.380 ±0.006 




„R 
BB 




0.469 ± 0.004 


0.532 ±0.006 


Mqq(m ICC ) 


f A 

c 

/"comb 
^comb 


Argus fraction 

Argus cndpoint (GeV/c 2 ) 

Argus exponent 

erf fraction 

erf width (GeV/c 2 ) 


1.0 

1.8696 
-17±9 
0.52 ± 0.15 
0.0005 ± 0.0002 


1.0 
1.8696 
-17 
0.52 
0.0005 


■Mcomb(m re c) 


f Al 

Cl 
C2 
tqq 

(Jqq 


Argus fraction 

A tcfuq ptiH Tiriint ( f-lp\/ / \ 

Ai exponent 
A2 exponent 
erf fraction 
erf width (GeV/c 2 ) 


0.27 ±0.06 
1.8696 
-49 ± 38 
-0.56 ±0.25 
0.26 ±0.04 
(75 ± 9) ■ 10 -5 


1.0 
1.8695 
-15 ±2 

0.41 ±0.06 
(72 ±2) ■ 10~ 5 


M B i g {m rcc ) 


fa 
m G 
OG 


Gaussian fraction 
Gaussian peak (GeV/c 2 ) 
Gaussian width (GeV/c 2 ) 


0.46 ±0.06 
1.8638 ± 0.0002 
0.0017 ±0.0002 


0.64 ±0.12 
1.8635 ± 0.0003 
0.0019 ± 0.0003 



^2500 




.6000 



1.84 1.845 1.85 1.855 1.86 



1.865 1.87 

m rec (GeV/c 2 ) 



FIG. 3: Result of the kinematical fit of kaon (top) and lepton 
(bottom) tagged data events, with PDF's overlaid: total PDF 
(solid line) , total background (BB + continuum, short-dashed 
line), BB combinatorial background (dashed line), continuum 
u, d, s, c background (dot-dashed line) and signal (dotted red 
line). 



To reduce dependency on Monte Carlo the fraction a of 
£>-tags remaining after the cos #tag constraint is measured 




g 1000- 



FIG. 4: Result of the kinematical fit of kaon (top) and lepton 
(bottom) tagged events for the Fisher function, with PDF's 
overlaid: total BB (solid line), BB (dashed line) and contin- 
uum it, d, s, c background (dash-dotted line). 



using data, as explained below. 

Defining the total number Nt (Nd) of T (D) tags, 
the number Nj, (N l D ) of them that satisfy the cos#t a g 
requirement, and the number Nip (Nfj) of them that do 
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not, a is given by: 

a = N l D /{N l D + N l T ) = 



f D P D 

(f D P D + (i-f D ) P T r 



(20) 



where p T = N l T /N T {p D = N l D /N D ) is the probability, 
taken from signal Monte Carlo, for a T-tag (D-t&g) to be 
from a track satisfying the cos6* tag cut. The fraction of 
D-tags, f D = N D /(N T + N D ), is given by 



f D = (p 1 -/W -p u ) 



(21) 



The fraction f = N l /{N l + N 9 ) is obtained from the 
kinematic fit of the data: N l is the number of signal 
events that have cos(9 tag < 0.75 (0.50) and N 9 is the 
number of signal events with cos6> ta g > 0.75 (0.50) for 
kaon (lepton) tag events. 

In this way we obtain the values a = 0.12 ± 0.04 for 
kaon tags and a — 0.00 ± 0.02 for lepton tags, as shown 
in Table |Hlj where we also list the mistag parameters oj 
and Aw, a, Tf,, and Amj that we will need in the final 
At fit. 

We use information from signal Monte Carlo events 
to determine the mistag probability oj = 0.201 ± 
0.002 (0.101 ± 0.002) and mistag difference Auj = 
-0.011±0.003 (0.001±0.005) for the kaon (lepton) tagged 
samples. We use the world average value for Amj [18], 
and the values fitted to signal Monte Carlo for Tf,, as 
discussed in Sec. lIIIBl 

TABLE III: Values of mistag parameters and a used in the 
final fit. The b lifetime values were obtained from the fit of 
signal Monte Carlo. The statistical uncertainties are given. 



parameter 


kaon tags 


lepton tags 


UJ 


0.201 ±0.002 


0.104 ±0.002 


Auj 


-0.011 ±0.003 


0.001 ± 0.005 


a 


0.12 ±0.04 


0.0 ±0.02 


Tb (PS) 


1.458 ± 0.014 


1.518 ±0.018 


Am d (ps _1 ) 


0.507 ±0.004 



the data strongly disfavor a component with nonzero life- 
time, therefore we fix fg = 1, and leave free in the final 
At fit only Aojs from Eq.[l5] For lepton tags, as the real 
data sample is too small, we obtain resolution and At 
parameters from continuum Monte Carlo. 

We use the continuum parameters obtained above in 
the next fit stage, where we obtain the BB background 
resolution function and At parameters using real data, by 
restricting the fit to events in a sideband region (1.836- 
1.856 GeV/c 2 ) of the D° recoil mass distribution. Ac- 
cording to Monte Carlo studies the fraction of signal 
events in this sideband is negligible, and we set it to 
zero. We fit simultaneously the resolution and lifetime 
parameters, shown in sections 1 and 3 of Table |IV| The 
fitted values of C com b and S^omb are consistent with as 
expected. 

We are now in a position to perform the final At fit, 
in which we release parameters from the signal, contin- 
uum and BB combinatorial background At and resolu- 
tion models, as detailed in Table |IV| always with the 
convention that parameters free (fixed) in the final At fit 
are shown with (without) a fit uncertainty. 

The results are also shown in Fig. [5] ^ for the kaon 
(lepton) tagged samples, where we plot the At distribu- 
tions separately for B° and B° tags, together with the 
time-dependent raw CP asymmetry 



A{At) 



N Sta&=1 (At) - Ns^.^At) 
N^-^At) + N St -^(Aty 



(22) 



For presentation purposes, only data in the restricted 
signal region m rcc > 1.860 GeV/c 2 are shown in Fig. [5] 
and Fig. [6] in order to amplify signal/background ratio 
and be able to see the oscillation. The signal fractions 
in this region become « 24% and « 18% for kaon and 
lepton tags respectively. 

This requirement is not applied to the fit sample, so our 
numeric results apply to the whole signal region m rcc > 
1.836 GeV/c 2 . 



SYSTEMATIC UNCERTAINTIES 



C. Time Dependent Fit 

After the kinematic fit is complete and mistag proba- 
bilities are determined, we can proceed to the final At fit 
to extract the CP-violating parameters S and C. 

We start by fitting the signal At PDF and its reso- 
lution function using a pure signal Monte Carlo sample; 
the parameters determined in this way will be used to 
initialize the signal PDF, and some of them will be left 
free again in the final At fit. 

Next we fit the resolution function parameters and 
the effective lifetime of the continuum background, us- 
ing the off-peak data sample. For the kaon tag sample, 



Our systematic uncertainties on the CP-violating pa- 
rameters S and C are summarized in Table IVl We discuss 
here the most significant ones. 

Most systematic uncertainties in Table [V] are due to 
imperfect knowledge of one single parameter fixed in the 
final At fit, having little or no correlation with uncer- 
tainties of other parameters. They have been treated by 
varying them by ±1<7 and repeating the final At fit leav- 
ing only the parameters S and C free to vary. 

Uncertainties in the first two lines have a different char- 
acter because they are due to parameter sets, in which 
correlations among parameters belonging to one set are 
non-trivial. Given the low signal-to-background ratio, 
correct modelling of the background shape and signal 
fraction in the kinematic fit is crucial, especially because 
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TABLE IV: Result of the final full fit. The values of fixed parameters are given without uncertainties. 



PDF 



BB resolution model 



continuum At 



parameter 

b" 
b° 
b w 

r 

s n 
s° 
s w 



description 

offset of narrow Gaussian 
offset of outlier Gaussian (ps) 
offset of wide Gaussian 
fraction of narrow Gaussian 
fraction of outlier Gaussian 
see eq. 

see eq. TB (ps) 



see eq. T7 

see eq 
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kaon tags 

-0.16 ±0.01 
0.0 

-1.0 ±0.2 
0.93 ±0.01 
0.008 ±0.001 
1.03 ±0.03 
8.0 
3.0 

-0.04 ± 0.02 



lepton tags 

-0.022 ± 0.014 
0.0 

-0.7 ±0.7 
0.977 ±0.004 
0.006 ±0.002 
1.02 ±0.02 
8.0 
5.6 

0.03 



BB At 



signal resolution model 



f 

J comb 

AWcomb 
Teomb 
Ccomb 
Scomb 

b" 
b w 

r 
r 

s n 
s w 
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see eq 
see eq 
see eq. Tl 

effective lifetime (ps) 
cosine coefficient 
sine coefficient 

offset of narrow Gaussian 
offset of wide Gaussian 
fraction of narrow Gaussian 
fraction of outlier Gaussian 
see eq. : 



see eq. T7 
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0.10 ±0.02 
0.04 ± 0.12 
-0.025 ±0.012 
1.318 ±0.023 
-0.022 ±0.024 
0.004 ± 0.014 

-0.35 ±0.09 
8±3 
0.992 ±0.007 
0.0 

1.13 ± 0.12 
2.6 



0.25 ±0.02 
-0.08 ±0.07 
0.012 ±0.021 
1.272 ±0.031 
-0.024 ± 0.041 
-0.023 ± 0.024 

-0.3 ±0.2 
1.0 ±0.1 
1.17±0.21 



signal At 


C 


±0.117±0.111 


±0.195 ±0.147 


S 


-0.417 ±0.159 


-0.210 ±0.200 



the M(m ICC ) and T(F) PDF's parameters are fixed in 
the final At fit. Consequently, we devised a procedure 
to evaluate the associated systematic uncertainties, that 
would also preserve the correlations among parameters 
belonging to a set. 

For each set of parameters in A4 or T that become 
fixed at any stage of our fits, and are not released again in 
the final At fit, a large number N t of toy Monte Carlo ex- 
periments of the same size as the data are generated and 
fitted, and the values of parameters in the N t experiments 
are saved. Evaluation of the systematics due to a set 
of parameters in subsequent fits (in which they become 
fixed) is made by repeating the latter fits many times 
over, using the same event sample, but fixing parameters 
in the set to different values every time, taken from one 
of the Nt experiments. In this way, we propagate the 
variation associated to parameter sets from one fit to the 
next one, and preserve correct correlations among them. 
We applied this procedure to obtain the uncertainties in 
lines 1 and 2 of Table [Vj which for lepton tagged events 
are the main source of systematics. 

For lepton tags we find that only one Gaussian is suf- 
ficient to describe the resolution function (/" = 1). The 
systematic due to the signal s w was evaluated only for 
kaon tagged events. 

Since the mistag parameter to is obtained from Monte 
Carlo with a very small statistical uncertainty (see Ta- 



ble III), we verified the agreement between Monte Carlo 



and data using a control sample of self-tagging B° 
D* =F 7r ± events. As a result of this study, we assign a 
very conservative uncertainty of 15% on w, and evaluate 
the associated systematic by repeating the final At fit 
varying its central value by ±15%. This is the largest 
systematic uncertainty for kaon tagged events. 

We estimate the systematic uncertainty associated 
with fixing the peaking background fraction fa in Eq.[8] 
to zero by setting it to ±0.002 for both the kaon and lep- 
ton tag samples, repeating the fit, and taking the largest 
deviation from the value fitted with fa — as the sys- 
tematic uncertainty. 

The signal m roc spectra for the CP-even and CP-odd 
components arc different, the latter being slightly harder. 
This may cause a small acceptance difference of our event 
reconstruction and selection, leading to a systematic shift 
in the C and S measurement. We have carefully evalu- 
ated this effect and found it to be negligible. 

As the D* + D*~ final state is a superposition of CP- 
even and CP-odd wave functions, the measured values S 
and C from our data only represent a weighted average 
of these components, with their inverse squared errors 
as weights. Since the background shape is not uniform 
as it goes to zero at the kinematical limit, the weight 
of the CP-odd component could be enhanced with re- 
spect to the CP-even one by the lower background level 
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800 



At (ps) 



FIG. 5: Top: At distribution for B° (dashed) and B° (solid) 
kaon tags; the lower curves are the corresponding signal 
PDFs. Bottom: raw time-dependent CP asymmetry. Only 
data in the restricted signal region m r . 
shown. 




> 1.860 GeV/c 2 are 



At (ps) 



FIG. 6: Top: At distribution for B° (dashed) and B° (solid) 
lepton tags; the lower curves are the corresponding signal 
PDFs. Bottom: raw time-dependent CP asymmetry. Only 
data in the restricted signal region m roc > 1.860 GeV/c 2 are 
shown. 



in the high mass region. To evaluate this effect, we per- 
form At fits in the two extreme Monte Carlo configura- 
tions, adding to the background sample a pure CP-odd 
(R± = 1) or CP-even (R± = 0) sample of signal events, 
respectively. The number of signal events in both cases 
is equal to the number of signal events found in data. 
We find that the differences in the errors of S and C 
are negligible in these two cases and we do not assign a 
systematic uncertainty to this effect. 

As discussed in Sec. |IV C"l Ti and Arrid are fixed to 
the values listed in Table |Hlj We assign the systematic 
uncertainty due to these assumptions by varying their 
nominal values of ±lcr, and taking half the difference in 
the fitted values of C and S so obtained. 

To evaluate bias on C and S in our fit, we apply the 
fit procedure to pure signal Monte Carlo events and com- 
pare the results for C and S to the generated ones; no 
significant bias has been observed in either. We therefore 
quote the statistical uncertainty on these Monte Carlo 
measurements as the associated systematic uncertainties. 

To measure the systematic uncertainty related to im- 
perfect knowledge of the time measurement due to uncer- 
tainty in the boost or possible uncorrected misalignment 
of the SVT, we repeat the time-dependent fit with dif- 
ferent sets of realistic misalignments of the SVT and At 



scaling factors. The maximum observed shift with re- 
spect to the nominal fit is quoted as the uncertainty. 

An important source of systematic uncertainty in our 
analysis is represented by interference effects from dou- 
bly Cabibbo-suppressed decay amplitudes on the tagging 
side of the event. The non-leptonic B-meson decays used 
for tagging are dominated by amplitudes containing a 
b — > cud transition, thus ensuring the correlation of the 
tagging particle (typically a kaon) with the flavor of the 
originating b quark. However, b — > Tied transitions could 
also contribute, although they are suppressed [2 2] by a 
factor r> * \{V* b V cd ) / {V cb V* d ) \ = 0.02. 

As discussed in detail in Ref . [35] , this effect cannot be 
simply reabsorbed into the mistag probability uj because 
the allowed and doubly Cabibbo-suppressed amplitudes 
can interfere, and thus effectively change the At proba- 
bility density function. 

Since our At PDF assumes r' = and therefore does 
not include these effects, the C, S parameters measured 
by our fit will be different from the observables without 
tag-side interference by a calculable quantity. 

To evaluate the systematic effect in our measurement 
due to neglecting small terms in the PDF with r' ^ 0, we 
follow the prescription in Ref. |22) and perform a simple 
toy Monte Carlo of SC = C fit - C and 8S = S m - S , 
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TABLE V: Systematic uncertainties evaluated for C and S. 
Uncertainties in the top section are independent for kaon and 
lepton tags, those in the bottom section are correlated. 





kaon 


tags 


lepton tags 


Systematic source 


C 


S 


c 


s 


Kinematic fit parameters 


0.013 


0.034 


u.Uzo 


U.UD / 


Continuum At fit parameters 


0.002 


0.001 






Signal s w 


0.0002 


0.0007 






BB combinatorial s w 


0.017 


0.0007 


0.001 


0.005 


Signal tag side (u>) 


0.012 


0.045 


0.002 


0.002 


Mistag difference (Acj) 


0.007 


0.0004 


0.007 


0.0009 


Signal CP side (a D o) 


0.006 


0.017 


0.002 


0.002 


Peaking background 


0.0002 


0.0003 


0.0002 


0.00004 


Fit bias (MC statistics) 


0.011 


0.018 


0.012 


0.019 


Tag interference from DCSD 


0.030 


0.002 






B a lifetime variation 


0.0002 


0.002 


0.0003 


0.004 


Amd variation 


0.0003 


0.001 


0.0004 


0.002 


SVT misalignment 


0.003 


0.007 


0.002 


0.004 


Boost uncertainty 


0.002 


0.006 


0.005 


0.007 


Total 


0.042 


0.062 


0.028 


0.061 



and CP odd wave function components. If penguin am- 
plitudes can be neglected then S + = —S-, C + = — C_ 
and the value of the CP-even components S + and C + , 
which we are interested in, can be obtained using the 
relations: 

C = C+ 

S = S+(1-2R X ), 

where the factor (1 — 2R±) represents the dilution in- 
troduced by the CP-odd component R± in the signal. 
To compute S + we use the value measured by BABAR of 
(R± = 0.158 ± 0.029) [5], where the uncertainty is the 
combined statistical and systematic. To evaluate the re- 
lated systematic uncertainty, we vary this value by ±lo\ 
We obtain 

C+ = +0.15 ±0.09 ±0.04 

S+ = -0.49 ±0.18 ±0.07 ±0.04, 

where the uncertainties shown are statistical and system- 
atic; the third uncertainty is the contribution from the 
error on R± described above. 



finding the results reported in Table |V| The lepton tags 
are not affected by this issue. 



VI. PHYSICS RESULTS 

The final results for C and S, with their correlation 
coefficient p, including only the statistical uncertainty 
for kaon and lepton tags, are: 



C = +0.12 ±0.11 
S = -0.42 + 0.16 



p = 0.0601, kaon tags, 



VII. SUMMARY 

We have measured the time-dependent CP asymme- 
try parameters C and S in B° — s- D* + D*~ decays, from 
which we have extracted the CP-even components S+ 
and C+. This result is an independent determination 
of the CP- violating parameters of b —> (cc)d transi- 
tions and is compatible with previous measurements from 
BABAR [5] and Belle [S] using fully reconstructed decays. 
It also agrees well with the Standard Model expectation 
of negligible contributions to the decay amplitude from 
penguin diagrams and thence with S + = — sin 2(3. 



C =+0.20 + 0.15 
S = -0.21 + 0.20 



p = 0.0730, lepton tags. 



The two samples are statistically independent of each 
other and can therefore be combined; their statistical 
uncertainties can be combined in quadrature, but the 
systematic ones need a more careful treatment. 

Indeed, several of the systematic effects listed in Ta- 
ble|V] are independent for the kaons and lepton tags and 
are combined in quadrature, while the others are com- 
bined taking into account their correlation. Finally we 
get the combined results of this analysis of 



C = +0.15 + 0.09 + 0.04 
S = -0.34 + 0.12 + 0.05 



0.0649. 



A. Extraction of S+ and C+ 

The measured values of S and C that we obtain from 
data only represent a weighted average of the CP even 
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